Gauge Dependence of the Effective Average 
Action in Einstein Gravity 



Sven Falkenberg 

Universitdt Leipzig, Institut fur Theoretische Physik 
Augustusplatz 10, D-04109 Leipzig, Germany 
e-mail: Sven.Falkenberg@itp.uni-leipzig.de 

Sergei D. Odintsov 

Dept. de Fisica, Universidad del Valle 
A. A. 25360 Call, Colombia 
and 

Dept. of Physics, Hiroshima University 
Higashi- Hiroshima, Japan 
on leave from State Tomsk Ped. Univ., 634041 Tomsk, Russia) 
e-mail: odintsov ©quantum, univalle. edu. co 
sergei@ecm.ub.es 



Abstract 

We study the gauge dependence of the effective average action and 
Newtonian gravitational constant using the RG equation for IV Then 
we truncate the space of action functionals to get a solution of this equa- 
tion. We solve the truncated evolution equation for the Einstein gravity 
in the De Sitter background for a general gauge parameter a and obtain 
a system of equations for the cosmological and the Newtonian constants. 
Analysing the running of the gravitational constant we find that the New- 
tonian constant depends strongly on the gauge parameter. This leads to 
the appearance of antiscreening and screening behaviour of the quantum 
gravity. The resolution of the gauge dependence problem is suggested. For 
physical gauges like the Landau-De Witt gauge the Newtonian constant 
shows an antiscreening. 



1 Introduction 



Despite many efforts during the last decades a consistent theory of quantum grav- 
ity (QG) is still lacking. The standard Einstein gravity is known to be non- 
renormalizable. But also renormalizable models like i? 2 -gravity (see jl] for a re- 
view) have other serious problems instead of this (e. g. presumably non-unitarity 
of the S-matrix). 

In such circumstances it is quite reasonable to consider some of the known 
pre-QG theories as effective models which appear from fundamental theories and 
which are valid below a definite scale only (e. g. Planck mass). In such a spirit 
one can introduce models to describe QG in the far infrared (at large distances) 
0, |3|, |J, new models of four- dimensional QG like dilatonic QG [|] and so on. 
For example, taking account of IR effects in QG with cosmological constant an 
effective model for the conformal factor has been constructed. Due to strong 
renormalization effects in the infrared this model may provide the solution of the 
cosmological constant problem. On the other hand, the type of theories to de- 
scribe the effective IR dynamics in QG maybe extended to quantum i? 2 -gravity 
DU and to quantum supergravity |§. Due to the fact that the effective conformal 
factor model maybe formulated as renormalizable theory one can calculate 
the (perturbative) running of the Newtonian coupling constant in such a theory 
0. Note that the running of cosmological and Newtonian coupling constants 
also can be discussed in other renormalizable theories (GUTs in curved space- 
time, i? 2 -gravity ||) where results of one- loop calculations of the corresponding 
/3-functions are available. Of course, all such considerations are closely related 
with the perturbative approach to the effective action in QG. 

However, if one still insists in using the non-renormalizable Einstein theory 
as low-energy effective theory for QG other methods should be employed. In 
particular, the use of non-perturbative RG in Wilson's effective action approach 
II maybe quite useful. In such an approach the non-renormalizability of the 
theory is no problem since the RG equation governing the evolution of the effective 
action from a scale A, where the theory is well defined, to a smaller scales k < A 
can be constructed. Such a finite RG equation (or evolution equation) maybe 
consistently formulated in the frame of the effective average action (T^, [II], [12], [U| 
(for recent discussions on related non-perturbative RG approaches see O] and 
references therein). 

In a recent work JT5| the scale-dependent gravitational average effective action 
Tfcf^y] has been discussed in the background field formalism. The corresponding 
exact and truncated evolution equations have been derived, and a special solution 
of the truncated equation has been given. The scale dependence of the cosmolog- 
ical and Newtonian constants in Einstein gravity has been determined (for earlier 
discussions of the average effective potential in QG see JTB[] ) . The whole study has 
been done using the background harmonic gauge (Fock-De Donder gauge). It has 
been found that the Newtonian coupling constant increases at large distances. 

On the same time it is quite well-known that the standard effective action in 
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the background field method depends on the choice of the gauge fixing condition 
(for a review, see |l[]). This gauge dependence is controlled by the corresponding 
background Ward identities (Slavnov- Taylor identities). Moreover it is usually 
assumed in quantum field theory that the effective action (if non-perturbatively 
calculated) does not depend on the gauge condition at its stationary point. Of 
course, in practice where only one-loop results are easily available it often happens 
that the effective action even at stationary points depends on the gauge fixing. 
This could lead to an apparent (or intermediate) gauge dependence also of some 
physical quantities. 

Similarly, in the study of effective average action with the background field 
method its gauge dependence is obvious. Such a gauge dependence is controlled 
by the background Ward identities again [[H]]. However it means that results 
obtained by investigation of the evolution equation (especially in the truncated 
form) necessarily should depend on the gauge choice. Hence it is interesting 
to investigate this gauge dependence for a scale-dependent gravitational action 
Tk[9nv] and for the associated evolution equation. 

In the present paper we discuss this question of gauge dependence for Ein- 
stein gravity with cosmological constant in a class of gauges depending on some 
parameter a where a = 1 corresponds to the harmonic gauge. 

The paper is written as follows. In the next section, following [[15], we briefly 
review the scale-dependent gravitational action and derive the exact and truncated 
RG equations. Section 3 is devoted to the calculation of the one-loop partition 
function for Einstein QG on De Sitter background (we work in the a-parametric 
class of gauges). This one-loop partition function actually defines the truncated 
evolution equation. In section 4 the evolution equations for the cosmological 
and gravitational constants are explicitly obtained with account of the gauge 
dependence. The running of the Newtonian coupling constant is analysed in 
section 5. It is shown that for some region of the gauge parameter Einstein 
gravity shows antiscreening (a > —2, 73) while for other values of a it shows a 
screening effect. The resolution of this unphysical gauge-depending situation is 
presented in the final section. 



2 Exact and Truncated Renormalization Group 
Equation 

We start with a short review of the average effective action formalism in QG; 
mainly we follow the exposition given in |15 . 

First we introduce a scale-dependent modification of the generating functional 
for connected Green's functions of QG in four dimensions 
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Jvh^Va i VC^ex V {-S[g + h] - S gi [h;g] (2.1) 

—S g h[h, C,C;g] — A k S[h, C,C;g] — Source} . 

Here we decomposed the integration variable 7^(2) being the gravitational field 
into a fixed background metric g^ u and into h^ u , the small fluctuations around 
this background, 

l^v{x) = g^u(x) + h^ v (x) . (2.2) 

This background field method (for an introduction, see permits to replace 
the integration over 7 by an integration over h. and are the ghost and 
anti-ghost fields for the gravitational field, respectively. 

The first term of the action, S[g + h], denotes the classical part being invariant 
under general coordinate transformations £7 generated by the Lie derivative with 
respect to an arbitrary vector field. At the moment we assume that S is positive 
definite, since the functional integral makes sense only if it converges. However, 
as has been argued in [15| this assumption is not necessary for the validity of the 
evolution equation which we want to derive (also the use of Lorentzian metric 
changes this equation only in an obvious manner). 

The gauge fixing term Sgf is given by 

S gi [h;g] = ^ fd'x^r v F,F u (2.3) 



with 



F, = V2KTf[g) h a(3 = V2k ^g^D, - \g af3 D^ K t 



We introduced here the constant k = (32ttG) 2 where G is the bare Newtonian 
constant. The derivatives are covariant with respect to the background metric 

By S g h we denote the ghost action which can be written for this choice of 
gauge fixing with the help of the Faddeev-Popov operator 

M [7, gfu = r p g aX Dx(l P »D* + lavDp) - g pa r X Dx 1* V D P 



S gh [h,C,C;g} = -V2jd i x^C fl M[g + h,gr u C ! ' . (2.4) 
That part of the action introduced so far, namely S + S g f + S g f , is invariant under 



BRST transformations which will not be specified here (see |T3j ) . 

In addition we introduced an IR cutoff for the gravitational field and for the 
ghosts: 
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A k S[h,C,C;g] = X -k 2 fd*x yfi h flv B** r \g]'»"»h pi , 

r (2.5) 

+V2 d 4 x^C^Rt h [g}C» . 



-Rf rav [g] and -Rf h [g] are the cutoff operators being defined below. They discriminate 
between low- and high-momentum modes. Whereas eigenmodes of —D 2 with 
eigenvalues p 2 ^> k 2 are integrated out in ( |2.1| ) without suppression, modes with 
eigenvalues p 2 <C k 2 are damped by a kind of momentum dependent mass term. 
This ensures the renormalizibility of the theory, which is an effective one at scale 
k. Both cutoff operators have the structure 

R k [g]=Z k k 2 R^(-D 2 /k 2 ) . (2.6) 

The dimensionless function R^ will be chosen as (see |I5|, [It], [TTJ, 



it interpolates smoothly between lim i^ ) (w) = 1 and lim R^(u) = 0. The 



factors Z k are different for the graviton and for the ghosts. For the ghosts Z k 



Z\ is a pure number but for the graviton Z k = Zf w is a tensor constructed from 

it is 

The last term in (p. II) is the source term 



the background metric. In the simplest case it is (2| rav )^' )cr = g pp g ua Zf :a ' v 



Source = - Jd^^irh^ + CJ^ + a^C, 

+P flu C c (g^ + V) + r^C v d v C^ 



(2.8) 



where Cc denotes the Lie derivative with respect to the ghost field. Obviously, 
the sources (3^ u and r M couple to the BRST variations of and C p . 

Now, to obtain the effective average action T^, we perform a Legendre trans- 
formation of W k . Therefore, let us first introduce the k -dependent classical fields 
(also called mean fields) 

Then the effective average action is obtained as 

r k [g,g,Z,t,l3,T] = !d A x^ {^h^ + o^ + 0%}-W h [t,wP,T-,g\ 

J _ (2-10) 

-A k S[g - g, £, f; g] . 

The first two terms correspond to the usual Legendre transform. By definition 
the effective average action is the difference between the usual Legendre transform 
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and the cutoff action A k S with the classical fields inserted. For k — > the cutoff 
A k S disappears and we get back the usual Legendre transform. In analogy to 7 M „ 
we introduced the metric 



9tiu(x) = g^{x) + hn U (x) 



(2.11) 



and we considered as functional of g instead of h. The effective average action 
T k depends on the mean fields as well as parametrically on the sources (3 and 
r and on the background metric g. It should be noted that this functional by 
virtue of the background field technique is invariant under general coordinate 
transformations. In this effective average action excitations with momenta larger 
than k are integrated out, while those with momenta smaller than k are suppressed 
by AkS such that T k describes an effective theory at scale k. In this sense the 
flow of Tfc with respect to k interpolates between an action containing the classical 
action T^oo = S+S g f +S g h — Jd^x -Jg {(3^ C^g ^ v +t^ u d v ^} = S and the common 
effective action = T. 

The trajectory in the space of all action functionals can be obtained as the 
solution of renormalization group equations. S enters via the initial condition 
for the trajectory at some UV cutoff scale A: T\ = S. While for renormalizable 
theories A may be send to infinity at the end, for QG it is a fixed, finite scale. 

Having now introduced the effective average action we derive the evolution 
equation for T k according to the procedure introduced in [II], |TTl [T2|, [13], [T5[ . 
We take the derivative of the eq. (|2.1|) with respect to the renormalization group 
"time" t = Ink and, after Legendre transformation, express the RHS in terms of 
T k using the relation between connected and 1PI two-point functions. Thus we 
obtain the following exact evolution equation 



d t T k 



-Tr 



- -Tr 

2 



+ ) _ _ ( d t Rk 

/ hh V J hh 



lf> + R 



/ ) — ( r[ 2 ' ) + Rk 



dtRk 



(2.12) 



(2) 

where T k denotes the second functional derivative of T k , e.g. 



,(2) 



5T, 



(2.13) 



Rk is a matrix in the field space whose non-zero entries are 



\ fiupa 



Rk)__ 

/ hh 

Rk 



<S(R, 



_ = V2Rt[g] = -(R k 



(2.14) 



The traces in ( |2.12j) have to be interpreted as Jd 4 x yfg ; they are convergent, 
in the UV- and in the IR-region, for all cutoffs which are similar to (|2.7|) . The 
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main contributions come from the region of momenta around k, because of the 
appearance of the factor d t R k ; large momenta are exponentially suppressed. 

Now the problem is to solve eq. fl2.12| ). One way to find at least approximate 
but still nonperturbative solutions is to truncate the space of action functionals. 
This means that we work on a finite-dimensional subspace parametrised by only 
a few couplings. Let us in a first step try to neglect the evolution of ghosts. This 
leads us to the following ansatz 



r*[0,0,£,£;/8,T] = r k [g)+F k [g,g] + S s£ [g-g;g)+Sg l [g-g,£,£;g] 



Jd 4 x^{r^g, u + r,Cd v e} ■ 



(2.15) 



The classical ghost and gauge fixing terms have the same functional form as for 
the bare action. Also the couplings to the BRST variations do not change. The 
other terms depend on g and g. The first term 



f k [g] =r*[0, 0,0,0; 0,0] 



(2.16) 



is the most interesting one, while f k contains only the deviations for g ^ g such 
that V k [g, g] = 0. This leads to the interpretation of T k as describing the quantum 
corrections to the gauge fixing term Sg which also vanishes for g = g. 

It is easy to show that for an UV cutoff scale A — * oo the action fulfils 
the initial condition T\\g] = S[g] with S being the classical action. If we solve 
eq. ( P-12|) directly and analyse the initial condition for we find that I\ = 
must be satisfied. 

This suggests to set 1^ = for all k in a second step. From the Slavnov- Taylor 
identity of T k it follows (for details see that T k contains higher loop effects 
such that this kind of approximation really makes sense. 

Inserting ansatz ( |2.15| ) into the evolution equation (2.12) we obtain the trun- 
cated evolution equation which we wanted to derive 



d t T k [g,g] 



2 



-Tr 



K- 2 r k 2) i g ,g}+Rrw) d t Rrm 

-M[g,g]+R e k h [g]Y 1 d t Rt h [g} 



(2.17) 



Here we introduced 

T k [g, g] = V k [g\ + S gi [g -g;g\~ T k [g, g, 0, 0; 0, 0] 



(2.18) 



while r^ 2 -* is the Hessian of T k [g,g] with respect to g^ v at a fixed background 
metric g^ v . As a result of our ansatz the LHS of (|2.17|) does not contain any 
contributions from the ghosts. 
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3 One-Loop Quantum Gravity on the De Sitter 
Background 

In this section we will find the one-loop partition function for Einstein QG. The 
classical action is written as 



(3.1) 



where non-renormalizable Einstein gravity is considered to be valid below a UV 
scale A. The truncation 



G — > Gh 



G 



J Nk 



a 



a 



J Nk 



(3.2) 



is supposed to be done. In general also a depends on k. In this case we would 
have three /c-depending parameters instead of two. This would lead to a more 
complicated system of equations for these parameters which could be solved with 
numerical methods only. For the pure Yang-Mills theory such a calculation has 
been done We restrict ourselves to the case that a is a constant. 
Substituting ((Op in ( gTg ) we easily get |T§ 



(3.3) 



T k [g,g] = 2K 2 Z Nk Jd A x^ [-R(g) + 2X k ] 

+ ^z Nk Jd 4 x ^ r v (rfgap) (^r<v) • 

For the RHS of the evolution equation we need the second functional derivative 
of T k [g,g}. In order to calculate it we should specify some curved background. A 
useful choice is the De Sitter space which is maximally symmetric. This means 
that Riemann and Ricci tensors can be written as 



R 



pupa 

Rpu 



\dpp9ua 9ficr9up\R 5 
-;9puR , 



(3.4) 



where R is the curvature which works as an external parameter characterising the 
space. Using the decomposition 



hf.u = + - g^ v ip , g^ u = 



we find 



r (2) 

k grav 



2Z Nk k 2 jd A x \fg { - h,,,, 
1 



_□ _ 2A fc + -R 



-□ - 2X k 



D u h 



(3.5) 



flU 



(3.6) 
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,(2) 
fcrf 



2Z Nk K 2 — d 4 Xy/g 
a J 



D u h 



)J,V 



(3.7) 



However, only for the minimal choice a — 1, the non- minimal terms in the sum 
of (|3.6| ) and (|3.7|) cancel each other. For other choices of a we are dealing with 
non-minimal operators. In order to make a heat kernel expansion with respect to 
the curvature we need to work with canonical operators like A = □ + m 2 only. 
We solve this problem making a further decomposition of the quantum fields. 



Like in 19 we write 



(3. 



where = and D^h^ u = 0. It is necessary then to define the differentially 

constrained operators (if) is an arbitrary scalar) 



A (X)^ 



^ p (X)h 



[IV 



(-n + x)ip , 

g, v (-D + X) C u± 
Sffi (-□ + X) 



and 



(3.9) 



[IV 



Note that a similar technique has been discussed in [^0[ . 

Now the one-loop partition function in the gravitational sector can be easily 
calculated % [19], 



J grav 



Vh^Vip exp 



_ r (2) _ p(2) ' 
k grav fc gf 



det Ai( ~R(2a - 1) - 2aA fc 



det A ( -72 (a - 1) - 2a\ k 



2 



det A 2 ( -i2-2A jfc 



(3.10) 



detA (-2A fe ) 



Here, in the course of the calculations the decomposition ( |3.8|) is done and the 
corresponding Jacobian is taken into account (for more details, see [|19|j). 

As we see the one-loop partition function is expressed in terms of constrained 
operators. However, the transition to constrained operators according to ( |3.8| ) 
and (|3.9| ) introduces additional zero modes. This leads to a wrong answer when 
we expand the determinants of the constrained operators in fl3.10| ) as powers of 



the curvature using the heat kernel expansion. Therefore we will present Z grav in 
terms of unconstrained operators. It could be done using the following relations 
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det A V (X) = det(-D + X) v 
det A T (X) = det(-Q + X) T 



Finally, we obtain 



det Ai(X)det A [X 



■dot A (l I X-Ir 



R 
1 

dci Ao(A")(let A,.( .V - Ajj 



(3.11) 



v grav 









2 


det A T 


>- 


2A fe J 










1 

2 


det Ay 















det Ay ( - (2a - 1)R- 2a\ k 



detA 5 (-2A fc ) 



(3.12) 



In order to calculate the complete partition function we should also take into 
account the ghost contribution. The ghost action ( |2.4j ) simplifies for the De Sitter 
space and forg = g (this corresponds to a special choice of the metric field) to 
the standard form 



-□ - -R 
4 



(3.13) 



The corresponding partition function can be easily found: 



1 



det Ay — R 



(2) 



,(2) 



fcgrav> -fegf and the S host term 



(3.14) 



The complete one-loop partition function for the Einstein QG in the gauge 
is the product of ( |3.12j ) and fl3.14p where we set g — g. 

For the truncated evolution equation (|2.17|) we have to include the cutoff term 
AfcS" into our calculation. This means that we have to calculate = In where 
the path integral of Z^ contains in addition to V 
also A k S. 

The coefficients i?f rav and Zf h have to be chosen such that the kinetic term and 
the cutoff term combine to —D 2 + k 2 R^(—D 2 /k 2 ) for every degree of freedom. 
RS°) has to be seen as a dimensionless function. Because we did not make any 
decomposition of the ghosts and neglected all renormalization effects of the ghosts 
we can set Zf h = 1. Through the decomposition of the gravitational field into four 
parts the kinetic term of every part obtained an additional factor. These factors 
have to be compensated by Zf av containing projectors on every part. Also ^| rav 
is proportional to the renormalization factor Z^k- 
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Then the truncated evolution equation (|2.17[ ) reads (g = g) 

-i" 



d t Fk[g,g) = 2 Tr r 



with 



- Tr^ 



Af ( -□ + -i? - 2A fc + A; 2 i? (0) 



A/" ( -□ + ^-^i? - 2aX k + k 2 R^ 



-1 



AT f-0 + ^R- 2X k + A; 2 i? (0) ^ 



N \-n-2\ k + k 2 R ( - 0) 



-1 



'Nk 



M {z) = d t [k 2 R^\< 



Here the variable z replaces —D 2 /k 2 . 



(3.15) 



(3.16) 



4 Evolution Equations for Newtonian and Cos- 
mological Constant 



In this section we write down the renormalization group equation (|2.17|) for the 
action ( |3.3| ). We differentiate ( |3.3| ) with respect to t, then we set g^ v = g^ u . This 
corresponds to a special choice of a field theory. It follows that the gauge fixing 
term S g { vanishes such that the LHS of the renormalization group equation reads 

d t T k [g,g] = 2K 2 Jd 4 x^ [-R(g) d t Z Nk + 2 d t (Z Nk X k )] . (4.1) 

Here the dependence on the gauge fixing parameter a disappears. 

Now we want to find the RHS of the evolution equation. We expand the 
operators in ( |3.15| ) with respect to the curvature R because we are only interested 
in terms of order jd^x yfg and jd^x ^Jg R: 



A; 1 (bR-2a\ k + k 2 R^) = A: 1 (-2a\ k + k 2 R (0) ) 



(4.2) 



- Ar 2 (-2a\ k + k 2 R^)bR + 0{R 2 ) 
where a and b are free constants. For a more compact notation let us introduce 

A ia (z) = A,(-2aA fc + k 2 R^\z)) . (4.3) 
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These steps lead then to 



dtT k [g,g] 



MA T \ 



Tr 



v 



--Tr x 



2a - 1 



Tr, 



(4.4) 



+ 



KA V 2 \+0(R 2 ) 



The terms with JVq are the contributions of the ghosts. 

As a next step we evaluate the traces. We use the heat kernel expansion which 
for an arbitrary function of the covariant Laplacian W(D 2 ) reads 



Ti^Wi-D 2 )] = (4vrr 2 tr J (/)|g 2 [^] jd^Xy/g 

+ \Qx[W] Jd i x^R + 0{R 2 ) 



(4.5) 



by I we denote the unit matrix in the space of fields on which D 2 acts. Therefore 
trj(J) simply counts the number of independent degrees of freedom of the field of 
sort j = T, V, S, namely 

tr s (/) = l, tr y (/)=4, tr T (/) = 9 . 



The sort j of fields enters (|4.5|) via tr only. Therefore we will drop the index 
j of Aj a after the evaluation of the traces in the heat kernel expansion. 
The functionals Q n are the Mellin transforms of W, 



Qo[W] 

Qn[W] 



W{0) 

1 

Tin) 



dzz n ~ x W{z) , (n>0) 



(4.6) 



Now we have to perform the heat kernel expansion ( |4.5| ) in eq. (JOI). This leads 
to a polynomial in R which is the RHS of the evolution equation. By comparison 
of coefficients with the LHS of the evolution equation ( fOj) we obtain in order of 
Jd 4 x^g 



dti > ZNkXk) -ibjb) 



6Q 2 



NA~[ 



-i 



4Q 2 



MA 



-i 



8Q 2 



A^oAo 1 



(4.7) 
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and in order of fd^x Jg R 
1 1 



d t Z 



t^Nk 



24k 2 (4 



7T 



-30Q 2 
-12Q 2 
+ 4Qj 



"at a 2 

ATA- 1 



6 (2a - 1) Q 2 



ATA- 2 



AAAi 



-i 



(4.8) 



8Qi ATqAq 



-i 



We introduce the cutoff dependent integrals 

1 



r(n) y 
1 



~ dzz „-l« (0 K*)-*a (0), (* 



r(n) 7 



z 



n-1 



[z + R(°){z) +w\ 
[z + R(°)(z)+wf 



V 



(4.9) 



for n > 0. It follows that $o(w) = ®o( w ) ~ 0- + w ) P for n = 0. In addition we 
use the fact that 



J\f 



d t [Z Nk k 2 R^(-D 2 / k 2 )] 



^Nk 



[2 - r] N (k)} k 2 R^ + 2D 2 RW 



(4.10) 



with T)N(k) = — 9j(ln Z^k) = —dt^Z^k) / Z^ k being the anomalous dimension of 
the operator ^[g~R. Then we can rewrite the equations ( |4.7|) and (|4.8|) in terms of 
$ and $. This leads to the following system of equations 



dt(Z N k\k) 



An 2 (4 



7T 



6^(-2Ai/fc 2 ) +4$i(-2«VA; 2 ) -8^(0) 



3& 2 (-2\ k /k 2 ) +2& 2 (-2a\ k /k 2 ) 



(4.11) 



12k 2 (4 



7T 



- 30 $ 2 (-2A fe /A; 2 ) - 6(2a - 1) $ 2 (-2aA fc /fc 2 ) 
+ 6 $J (-2 A fc /A; 2 ) + 4 $J (-2a\ k /k 2 ) 
-12$ 2 (0)-8$K0) (4.12) 
-r] N (k) \- l5& 2 (-2\ k /k 2 ) -3{2a-l)& 2 (-2a\ k /k 2 ) 



+ 3 $l (-2 VP) + 2 $i (-2aA fc /fc 2 ) 

Now we introduce the dimensionless, renormalised Newtonian and cosmologi- 
cal constants 

g k = k 2 G k = k 2 Z N lG, X k = k~ 2 X k . (4.13) 
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Here G k is the renormalised but dimensionful Newtonian constant at scale k. The 
evolution equation for g k reads then 



d t 9k = [2 + f] N {k)}g k . 
From ( [4.12 ) we find the anomalous dimension rj^ik) 

r) N (k) = g k Bi(X k ) + r] N (k)g k B 2 {\ k ) 
where we used the notation 



(4.14) 



(4.15) 



£i(A fc ) 
B 2 {\ k ) 

Solving QH5|) 



±- [3 $}(-2A fe ) + 2 $l(-2aA fc ) - 4 $J(0) 

07T 



-15 $1(-2A*) - 3 (2a - 1) <&l(-2a\ h 
1 

67T 



6*1(0)] 



3$i(-2A fc ) + 2 4>l(-2aA fc ) 



(4.16) 



-15^(-2A fe ) 



VN{k) 



3 {2a - 1) $|(-2aAji 



g k Bi(\k) 



(4.17) 



1 - 9kB 2 (\k) 

we see that the anomalous dimension rj^ is a nonperturbative quantity. From 
( 4.11 ) we obtain the evolution equation for the cosmological constant 

d t \ k = - [2 - m {k)\ \ k + i- ^6 <^(-2A fc ) + 4 $^(-2aA fc ) - 8 $*(0) 



2tt 
-?7jv(/c) 



3&(-2A fc ) + 2&(-2aA A 



(4.18) 



The equations ( 4.14 ) and ( |4.18|) together with ( 4.17 ) give the system of differential 
equations for the two fc-depending constants X k and g k . If the initial values of g\ 
and Aa are known it determines the running of these constants for k < A. This 
system contains effects of arbitrarily high orders in the perturbation theory. In 
the special case a = 1 we get back the result of Reuter in K 



5 The Running of the Newtonian constant 

We will calculate now the a-dependence of the running Newtonian constant. To 
do this we start from the evolution equation for the renormalised gravitational 
constant 

d t G k = r ]N {k)G k (5.1) 



with ?in taken from ( J4.17 ) and g k = k 2 G k = k 2 G + 0{G 2 ). Let us assume that the 



cosmological constant is much smaller than the IR cutoff scale, X k <ti k 2 . So we 
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can neglect all contribution of X k in T]^ to simplify the result. Then the evolution 
equation has the solution 



G k = G [l-aGk 2 ] 



to/a 



with 



a 



B 2 (0) 



3a + 1 



CO 



B 1 (0) 



7i 

36 



36(a + 3) 



7T Z 



(5.2) 



(5.3) 



6vr 2 

2 

We see that u;>0foro;>-|g— 3~ —2.73. This means the gravitation is "anti- 
screening" in such a way that G k decreases if k increases, i. e. with increasing 
distance also the Newtonian constant increases. This supports the naive picture 
that the mass of a particle receives additional positive contributions from the 
virtual particles surrounding it. 

In addition there is a second interval for a where the antiscreening changes 
to a screening effect. The coefficient a does not play any role in the change from 
antiscreening to screening. The significant influence of the running of G k comes 
from the sign of the numerator of the exponent in ( |5.2| ) which is to. The value 
a = |g — 3 is significant for the used truncation and zeroth order in Afc. 

This screening behaviour is an artifact of the truncation where we used a(k) = 
a = const. A better truncation would include a as fc-depending parameter as it 
is done in \IE\ for pure Yang-Mills theory. There it has been shown that a = 
is an attracting fixed point of the truncated evolution. This is to be expected 
for gravity too. From ( |4 . 1 2| ) we see that Z^ k ~ k 2 . As in Yang-Mills theory we 
expect that a/Z^ k is a constant and does not depend on k. Then also a should 
tend to zero quadratically for k — > 0. 

To get back the result in [[15] we expand ( |5.2| ) with respect to Gk 2 and obtain 



G k = G [1 - w G k 2 + O (G 2 k A )] 



(5.4) 



with us given by ( |5.3j ) . For a = 1 Reuters result is obtained. 

In figure p] the gravitational constant G k is shown as function of the gauge 
parameter a and the cutoff scale k. The Landau-De Witt gauge a = is drawn 
as a separate curve on the surface. Both regions — antiscreening and screening 
region — are to perceive. 



6 Discussion 

In order to solve the gauge dependence problem we suggest two ways. The first 
way is the use of a better truncation which includes the gauge parameter as a 
function of k. As argued in the previous section a = should be a stable fixed 
point of the theory. It would be also interesting to use more general gauge fixing 
conditions. 

Alternatively it is possible to work with a generalisation of the standard effec- 
tive action, the so-called gauge-fixing independent effective action (for an intro- 
duction and its different forms, see [^l|, |l|]). The gauge-fixing independent effective 
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Figure 1: The gravitational constant as function of k and a. 



action leads to the same ^-matrix as the standard one; but by construction it does 
not depend on the choice of the gauge condition (off-shell it differs from the stan- 
dard effective action). Moreover, it has been shown |22| that in Einstein gravity 
on De Sitter background the one-loop standard effective action in the Landau-De 
Witt gauge a = coincides with the one-loop gauge-fixing independent effective 
action pi |. 



Working with the truncated evolution equation which is very similar to the 
one-loop approach one may give arguments that the truncated scale-dependent 
gravitational effective action at a = coincides with the gauge-fixing independent 
oneQ. Then the physical value of u would be given by ( |5.3|) at a = again. 

As one sees the Einstein gravity looks as an antiscreening theory where the 
gravitational constant increases at large distances: 



G(r) = G{oo) 



Gh (G 



(6.1) 



where r = k~ l , (G/r 2 <C 1), and the correspondence with the Newtonian potential 
may be restored. 

Of course, the behaviour of the gravitational coupling constant depends very 
much on the type of the model one considers. For example it has a qualitatively 
different form for the effective theory of the conformal factor 0, [7| , for i? 2 -gravity 
H or for the four- dimensional dilaton gravity M. In the very recent calculation of 
u in |23| (see also references therein) it was also found that u is positive. Similarly, 
in a = 1 gauge || uj is found to be positive again. 

Hence, our study of gauge dependence of truncated evolution equation in Ein- 
stein gravity indicates that the gravitational coupling constant increases at large 

1 Note that the latter effective action is not constructed yet as complete scheme. 



16 



distances. 

There are different possibilities to extend the results of the present work. The 
next problem is the calculation using an improved truncation. One can also con- 
sider a more complicated theory like gauged supergravity on gravitational back- 
ground, to formulate there evolution equation for scale-dependent effective action, 
and to study the running of gravitational and cosmological constants. We again 
will find the problem of the gauge dependence of the evolution equation. It could 
be also of interest to discuss in details the cosmological applications of gravita- 
tional and cosmological constants running (see |2^|). 

Finally, one has to note that taking into account of combined gravity-matter ef- 
fects (unified in supermultiplet as in supergravities theories or just Einstein-matter 
type systems) may also change the behaviour of gravitational and cosmological 
couplings. 
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